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The scalar isoscalar form factor a{t) of the nucleon is calculated in the limit of a large number 
of colours in the framework of the chiral quark-soliton model. The calculation is free of adjustable 
parameters and based on an approximation justified by the small packing fraction of instantons in 
the QCD vacuum model, from which the chiral quark-soliton model was derived. The result for 
cr(t) reproduces all features of the form factor observed in previous exact numerical calculations in 
' the chiral quark-soliton model and in chiral perturbation theory, and agrees well with the available 

I phenomenological information. The Feynman-Hellmann theorem is used to study the pion mass 

, dependence of the nucleon mass and a good agreement with lattice QCD results is observed. The 

' use of the present method to extrapolate lattice data to the chiral limit is discussed. 
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. The pion-nuclcon sigma-tcrm ct^n is of fundamental importance for understanding chiral symmetry break- 

^ ' ing effects in the nucleon [1, 2]. It can be inferred from pion- nucleon scattering data [3] and the value is 

I found (T^N — (50— 70) MeV [4-9]. The pion-nucleon sigma-term has wide phenomenological impact in many 

CO ■ fields. E.g., it is related to partial restoration of chiral symmetry in the nuclear medium [10, 11] and to 

CO I central nuclear forces [12]. It is an important ingredient in the mass decomposition of the nucleon [13] and 

plays an important role in the searches for the Higgs boson [14] , dark matter and super-symmetric particles 
[15, 16]. Valuable insights into the cr^N-physics were provided from studies in chiral perturbation theory 
^ . [17-22], lattice QCD [23-25], and numerous chiral models, e.g. [26-40]. 

In this note the scalar isoscalar nucleon form factor a{t), which at zero momentum transfer t yields CTttN, 
is studied in the chiral quark-soliton model (xQSM) [41, 42]. The picture of baryons as chiral solitons of the 
pion field is justified in the limit of a large number of colours Nc [43]. Despite the fact that in nature Nc = 3 
Qh! does not seem to be large the xQSM has proven to describe successfully numerous properties of the nucleon 

O . [44-47] and to provide valuable insights. E.g., in this model the observation was made of the so-called 

Z)-term in generalized parton distribution functions [48] - a fundamental characteristics of the nucleon [49] 
which, when known, e.g., will provide information about the "distribution of strong forces" in the nucleon 
[50] . Another recent highlight made on the basis of the soliton picture was the accurate prediction of the 
exotic 0^ baryon [51] for which by now strong experimental indication is collected [52]. 
' The xQSM was derived from the instanton model of the QCD vacuum [54, 55] and incorporates chiral 

symmetry breaking. The field theoretical character of the model plays a crucial role in ensuring the theoretical 
consistency of the approach. Both (Tttn and <T{t) were already studied in the xQSM in [27, 28] and [29], where 
the model expressions were exactly evaluated using involved numerical methods. Instead the calculation of 
a{t) presented here is based on an approximation which is consistently justified by arguments from the 
instanton vacuum model. 

The approximate method not only considerably simplifies the calculation - allowing a lucid interpretation 
of the results and suggesting an appealing explanation why CTttN is so large. In addition it has methodical 
advantages over the exact calculations [27-29]. For in the xQSM a{t) is quadratically UV-divergent and 
as such strongly sensitive to the details of regularization. The virtue of the approximate method is that it 
yields a regularization scheme independent result. Another advantage is that the theoretical accuracy of the 
results - which in models often has to be concluded from the comparison to phenomenology - is here known 
in advance and under theoretical control. To within this accuracy results from the exact xQSM calculations 
[27-29] are reproduced, and a good agreement with phenomenological information and results from other 
approaches is observed without any adjustable parameters. 

Another advantage is that the approximation allows to analytically study the properties of a{t) in the 
chiral limit - which is a rare occasion in the xQSM. Finally, by exploring the Feynman-Hellmann theorem 
[56], the pion mass (m^r) dependence of the nucleon mass (Mn) is studied. An exact treatment of this issue 
in the xQSM would be numerically highly involved. The results for M-[<i{m,r) compare well to lattice QCD 



2 



results. This observation could be used to inspire ansatze for the chiral extrapolation of lattice data. 

The note is organized as follows. In Section II the form factor (j(t) is defined and briefly discussed. In 
Section III the model is introduced. The form factor is computed in Section IV and the numerical results 
are discussed in Section V. Section VI is devoted to the study of a{t) and the nucleon mass as functions of 
the pion mass. Section VII contains the conclusions. Technical details of the calculation can be found in the 
Appendices A and B. 



II. THE SCALAR ISOSCALAR NUCLEON FORM-FACTOR 

The nucleon sigma-term form factor a{t) is defined as the form factor of the double commutator of the 
strong interaction Hamiltonian with two axial isovector charges [1]. Disregarding a "double isospin violating 
term" proportional to (m„ — md){'4>u'4'u — i^di^d) and commonly presumed to be negligible the form factor 
can be expressed as 

a{t) w(p')u(p) = m (Ar(p')l (^Vi„(0)V«(0) + ^d(0)Vd(0)) |iV(p)) , t={p- p'f , (1) 

where m = \{mu + md). In Eq. (1) |Af(p)) is the spin-averaged state of the nucleon of momentum p 
normalized as (A^(p')|iV(p)) — 2p^5^^^{p — p'), and u(p) is the nucleon spinor with 'u(p)zt(p) = 2Mn. 
At zero momentum transfer the form factor yields the pion-nucleon sigma-term, i.e. 

a^N ^ a(0) = ^ (iV(p)l (^^.(0)V„(0) + ^d(0)Vd(0)) |iV(p)) . (2) 

The pion-nucleon sigma-term cr^rN is normalization scale invariant. The form factor a{t) describes the 
elastic scattering off the nucleon due to the exchange of an isoscalar spin-zero particle. It is not known 
experimentally except for its value at the (unphysical) Cheng-Dashen point t = 2m^. A low energy theorem 
[3] relates the value of a{t) at the so-called Cheng-Dashen point t — 2m^ to the (isoscalar even) pion nucleon 
scattering amplitude. The analysis of pion-nucleon scattering data yields 

(64 ± 8) MeV (1982) [4] 

nClr^-^ 1 ^ ^ (88 ± 15) MeV (1999) [6] , . 

'^K^rn^) -\ ± g) ^cV (1999) [7] 

(79 ±7) MeV (2001) [8]. 

The recent analyses tend to yield a larger value for (T(2m^) which can be explained by the more recent and 
accurate data [9]. The difference ct(2to^) — (t(0) has been calculated from a dispersion relation analysis [5] 

o-(2m^) - o-(O) = (15.2 ±0.4) MeV . (4) 

In Ref. [21] a similar result was obtained from a calculation in the chiral perturbation theory. From Eqs. (3, 4) 

one concludes 

cr„N =i (50 - 70) MeV . (5) 

The large value of cr^N has been and still is a puzzle [1, 2]. According to a "standard interpretation" ct^n 
can be related to the so-called strangeness content y = 2{N\^s'4^s\N) / {N\{il)u'tpu + '4>d'4'd)\N) of the nucleon 
as (1 — y)(TnN = ^, where a can be determined by means of chiral perturbation theory from baryon mass 
splittings, a- = (35 ± 5) MeV [17, 57]. The value in Eq. (5) then implies y ~ (0.3 - 0.4) while one would 
expect y ~ on the grounds of the OZI-rule. It is worthwhile mentioning that such a sizeable value of y 
naturally follows from a Goldstone-boson-pair mechanism [26]. 



III. THE CHIRAL QUARK SOLITON MODEL (xQSM) 



This section briefly introduces the notions required in the following. More complete presentations of the 
model can be found in, e.g., [27, 41, 42]. The xQSM is based on the effective chiral relativistic quantum 
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field theory of quarks, antiquarks and Goldstone bosons defined by the partition function [55, 58, 59] 

U = exp(iT"7r'') , t/^^ = exp(i75r'^7r'') = ^{U + U^) + ^{U - U'')-f5 ■ 



(6) 



In Eq. (6) M is the dynamical qnark mass, which is due to spontaneous breakdown of chiral symmetry and 
in general momentum dependent. U = exp(ir''7r") denotes the SU(2) chiral pion field and m the current 
quark mass, which explicitly breaks the chiral symmetry. In many applications m can be set to zero, but for 
certain quantities it is convenient or even necessary to consider finite m,. The effective theory (6) contains 
the Wess-Zumino term and the four-derivative Gasser-Leutwyler terms with correct coefficients [42]. It has 
been derived from the instanton model of the QCD vacuum [54, 55] and is valid at low energies below a scale 
set by the inverse of the average instanton size 



600 MeV . 
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In practical calculations it is convenient to take the momentum dependent quark mass constant, i.e. M{p) 
M(0) = 350 MeV. In this case is to be understood as the cutoff, at which quark momenta have to be 
cut off within some appropriate regularization scheme. 

It is important to remark that (Afpav)^ if' proportional to the parametrically small instanton packing 
fraction, i.e. with denoting the average distance between instantons in Euclidean space-time 



< 1 . 



(8) 



Numerically Pav/-Rav ~ 1/3. The parameterical smallness of this quantity played an important role in the 
derivation of the effective theory (6) from the instanton model of the QCD- vacuum [41, 54]. 

The xQSM is an application of the effective theory (6) to the description of baryons [41, 42]. The 
large- Af^c limit allows to solve the path integral over pion field configurations in Eq. (6) in the saddle-point 
approximation. In the leading order of the large- A^c limit the pion field is static, and one can determine the 
spectrum of the one-particle Hamiltonian of the effective theory (6) 



H\n) = En\n) , H = -ij'^j'^dk + -f^MW^' + j°m 



(9) 



The spectrum consists of an upper and a lower Dirac continuum, distorted by the pion field as compared to 
continua of the free Dirac-Hamiltonian 



Ho\no) = Er,,\no) , Ho = -h^l'^dk + J°M + , 



(10) 



and of a discrete bound state level of energy Eie^, if the pion field is strong enough. By occupying the 

discrete level and the states of the lower continuum each by Nc quarks in an anti-symmetric colour state, 
one obtains a state with unity baryon number. The soliton energy Esoi is a functional of the pion field 



E,„i[U] = nJ El,, 



{En — En 

E„<0 



(11) 



Esoi[U] is logarithmically divergent and has to be regularized appropriately, which is indicated in Eq. (11). 
Minimization of Es^ilU] determines the self-consistent solitonic pion field Uc- This procedure is performed 
for symmetry reasons in the so-called hedgehog ansatz 



7r"(x) = e° P(r) , C/(x) = cos P(r) + ir^e^ sin P{r) , 



(12) 



with the radial (soliton profile) function P(r) and r = |x|, = x/r. The nucleon mass Mn is given 
by Esoi[Uc\. The momentum and the spin and isospin quantum numbers of the baryon are described by 
considering zero modes of the soliton. Corrections in the 1/iVc-expansion can be included by considering time 
dependent pion field fluctuations around the solitonic solution. The xQSM provides a practical realization 
of the large- A^c picture of the nucleon [43]. 
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The self-consistent profile satisfies Pc{0) = — tt and decays in the chiral limit as at large r. For finite 

m it exhibits a Yukawa tail cx exp(— m„r)/r with the pion mass connected to m by the Gell-Mann- 
Oakes-Renner relation, see below Eq. (34). An excellent approximation to the self-consistent profile, which 
exhibits all those features (and is sufficient for our purposes) is given by the analytical "arctan-profile" 

P(r) = -2arctan(^:^(l + TO^r)e-'""''^ , i?soi = M'^ . (13) 

The quantity R^oi is referred to as the soliton size. It is related to the nucleon axial coupling constant 
gA = 1.25 as [42, 60] (note that the order of the limits cannot be inversed) 

lim ( lim r'Pir)) = -2Rl, = f • (14) 

The xQSM allows to evaluate without adjustable parameters nucleon matrix elements of QCD quark 
bilinear operators as 

(iV(p')|V^(0)rV;(0)|7V(p)) = CrM^Nc J2 f e'(p'-P)-$„(x)r$„(x) 1^^^+ . . . (15) 

n.occ 

where F is some Dirac- and flavour-matrix, Cr is a constant depending on F and the spin and flavour quantum 
numbers of the nucleon state \N) = |S'3,T3), and $n(x) = (x|n) are the coordinate space wave-functions 
of the single quark states \n) defined in (9). The sum in Eq. (15) goes over occupied levels n (i.e. n with 
En ^ Eigv), and vacuum subtraction is implied for En < i?iev analog to Eq. (11). The dots denote terms 
subleading in the l/A^c-expansion (which can be included but will not be considered in this work). The 
model expressions can contain UV-divergences which have to be regularized as indicated in (15). 

If in QCD the quantity on the left hand side of Eq. (15) is normalization scale dependent, the model result 
refer to a scale of 0{p~^), see Eq. (7). In the way sketched in (15) a large variety of static nucleon properties 
like form-factors, axial properties, etc., were computed (see [44, 45] for reviews). In Ref. [46] the approach 
was generalized to non-local quark bilinear operators on the left hand side of (15) which paved the way to 
the study of the quark and antiquark distribution functions [46, 47] and off-forward distribution functions 
[48]. The model results agree typically to within (10 — 30)% with experimental data or phenomenological 
information. 



IV. a(t) IN THE xQSM 

In this Section first the model expression for a{t) is discussed and the consistency of the approach is 
demonstrated. Next the UV-behaviour of a{t) is studied and the question of regularization is addressed. 
Thereby is defined and justified the approximation in which then a{t) is evaluated in the next Section V. 

Expression and consistency. The pion-nucleon sigma-term ct^n was studied in the framework of the 
xQSM in [27, 28] and the scalar isoscalar form-factor a{t) in [29]. In leading order of the large A^c limit the 
model expression for the form-factor a{t) reads [29] (in the SU(2) flavour sector) 

cj{t) = m N,Jd^xjoiV^t\^\) J2 *:W7°*n(x)L,g , (16) 



, occ 



where the Bessel function jo{z) = sinz/ z. In the large- iVc hmit the nucleon mass Mn = 0{Nc) while the 

nucleon momenta |p| and ]p'] = 0{N^) such that t = — (p' — p)^. Therefore (16) is vahd for \t\ M^. 

(Interestingly, in the case of electromagnetic form factors the model results agree well with data up to 

\t\ ~ 1 GeV^ [44].) Eq. (16) shows that a{t) — 0{Nc) in agreement with results from large-A^'c chiral 
perturbation theory [22]. 

In the xQSM one is in a position to derive the model expression for CTttN in three different ways: 

(TttN = lim (j{t) , (17) 

cr^N = m — — — , (18) 

CT^N = m / dx (e" +e'^ + e'' + e^){x) . (19) 
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The first method (17) consists in continuing analytically the form factor a{t) to t = 0. The second method 

(18) uses the Fcynman-Hellmann theorem [56]. The third method (19) uses the sum rule for the first moment 
of the flavour singlet twist-3 chirally odd distribution function e°'{x) [61].^ The three methods consistently 
yield 

a.N = m iV, /d^x ^ $;(x)70$„(x)|^^g . (20) 

n, occ 

The result in Eq. (20) immediately follows from the model expression (16) (recalling jo{z) 1 for z — > 
0). Relation (18) was used to numerically compute ct^n in [27], and explicitly demonstrated to yield the 
expression in Eq. (20) in [63]. The sum rule (19) was shown to be satisfied in the xQSM (and to yield (20) 
for CTttn) in [63, 64]. The fact that CTttn can consistently be computed in the xQSM in three different ways 
illustrates the theoretical consistency of the model. 

The UV-behaviour of a{t). In this paragraph the known result, cf. [27-29], will be rederived that the 
model expression for a{t) contains quadratic and logarithmic divergences, i.e. that it is of the form 

= «W + '°<Tf )... + {¥)' • 

where the coefficients ai{t) are UV-finite functions of t and A^^j is an UV-cutoff. A similar study was 
presented in Ref. [63], however, for the more involved case of the twist-3 distribution function e"(a;) which 
is related to CTttn by means of the sum rule (19). 

Let us separately consider the contributions from the discrete level and the negative continuum 

a{t) = a{t)uy + cT(t)cont , cj{t)k^ = m jd^x jo{V^t ]x|) ^l^{yi)^Hi^^{x) , (22) 

a{tUnt =mN, /d^x jo(v^|x|) *nW7°^n(x)|,,g . (23) 

E„<0 

For the discrete level the eigenvalue problem jf$iev(x) = £^iev^iev(x) can exactly be solved numerically, see 
e.g. [46], and with <l>iov(x) one obtains (j{t)\ev The contribution of the discrete level is always finite and for 
our purposes it is sufficient to note 

(7(t)iev = UV-finite. (24) 

The exact evaluation of the continuum contribution in (23) is far more involved. For that one cither can 
place the soliton in a finite 3-D box, discretize and make finite the spectrum of the free Hamiltonian (10) 
by imposing boundary conditions and diagonalize the Hamiltonian (9) in the basis of the free Hamiltonian 
states (Kahana-Ripka method [65]). Alternatively one can rewrite the continuum contribution (23) in terms 
of Green functions and evaluate those by means of phase shift methods (see, e.g., Ref. [66]). Both methods 
are numerically involved. 

Here we will use an approximate method - in Ref. [46] referred to as interpolation formula - which consists 
in expanding the continuum contribution (23) in gradients of the J7-field and retaining the leading order 
only. The interpolation formula yields exact results in three limiting cases |Vf/| ^ M, |Vt/| 3> M and 
|logi7| <C 1. Therefore one can expect that it yields useful estimates also for the general case. It was 
observed that this method approximates exact calculations in the model with good accuracy [46, 47]. 

Let us rewrite the continuum contribution in Eq. (23) (recalling the implicit vacuum subtraction) as 



, , /■ dw 



Sp 



jo(y^|x|)7°— T--(^-^o) 
uj + iH 



(25) 



^ Eq. (19) is correct in a formal mathematical sense [61]. However, the sum rule (19) is saturated by a <5-function at a; = in 
e"(a;) which means that the relation (19) is pra<;tically useless to extracted any information on ct^n from possible measurements 
of e" (x) in deeply inelastic scattering experiments, see [62] and references therein. 
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where the contour C is defined as going along the real w-axis and closed in the infinity in the upper half 

of the complex cj-plane. The original expression (23) is recovered by saturating the functional trace with 
the complete set of eigenfunctions of respectively H and Hq in (9, 10),^ performing the w-integration, and 
passing to the coordinate space representation. Expanding (25) in a series in gradients of the iZ-field 

oo 

<7(t)cont = V ^it)^^, , (26) 



k=0 



where the index k means that VU appears fc-times in o'{t)['^^^, wc obtain (see Appendix A) 

C a B{t) for fc = 0, 
'^(i)coL = i zero for fc = 1, (27) 

I UV-finite for fc > 2, 



with the constant a and the "form factor" B{t) defined as 

d^B SN^M 



(28) 



(27r)>| + M2 

B{t) = j d^x io(\/=i |x|) - itrF[/(x)^ . (29) 



The factor SNcMm is included into the definition of a for later convenience. The Euclidean integral in the 
constant a contains quadratic and logarithmic divergences. Combining the results in (24) and (27, 28) we 
see that the form-factor cr(i) has the UV-bchaviour (21) (and that a2{t) oc a\og{t) oc B{t)). 

Instanton motivated approximation. In non-renormalizable effective (low energy) theories the regulariza- 
tion procedure "keeps the memory" of the cutoff Acut- In such effective theories - for which (6) is an example 
- the cutoff has a physical meaning: It sets the scale below which the degrees of freedom of the effective 
theory may be considered as appropriate to describe the physical situation, and above which they may not 
be sufficient. In the effective theory (6) the cutoff 

Acut = 0(Pav ) • (30) 

Using (30) and the results of the previous paragraph we see that (j{t) can be written as, 

cj{t)=aB{t)-{l + 0{M^pl,)\ (31) 



with a and B{t) as defined in Eqs. (28, 29). 

Thus, in Eq. (21) the UV-finite contributions are parametrically strongly suppressed with respect to the 
UV-divergent terms by the instanton packing fraction due to Eq. (8). Since the xQSM was derived from the 
instanton vacuum model, it is consistent to use this argument based on Eqs. (8, 30) in this context. Eq. (31) 
defines the approximation in which ait) will finally be evaluated, after regularizing the divergent constant 
a in Eq. (28). 

Regularization. There arc several methods to regularize the divergent constant a in Eq. (28). Popular 
methods often used in exact model calculations are the Schwinger proper-time regularization (see, e.g., 
[44]) or the Pauli-Villars subtraction method (see, e.g., [63, 67]). The result will be sensitive to the chosen 
regularization method and, especially for a power divergent quantity, to the precise value of the cutoff - of 
which only the order of magnitude is known, cf. Eq. (30). 

However, the practical problem of how to regularize the constant a in Eq. (28) can be solved elegantly and 
in a regularization independent way. Since the quark vacuum condensate (V'V')vac is given in the effective 
theory (6) by [55] 



))|vac)= I 



(27r)4 p| + M2 



(32) 



I.e. Sp[. ..-{H-. Ho)] = Eall„("l • ■ • H- Eallno ("o| ■ • • |no) 
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the constant a. in Eq. (28) can be expressed as 

a = - m (V''!/')vac • (33) 

At first glance Eq. (33) is based on tlie mere observation that the same divergent integral appears in two 
different model expressions. However, as will be discussed in the next section, the relation (33) is not 
accidental from a physical point of view (and implies an interesting interpretation of ct^-n)- In the next step 
one can use the Gell-Mann-Oakes-Renner relation (with denoting the pion decay constant) 

mi!l = - m (V;V)vac , (34) 

which is not imposed here by hand but valid in the effective theory (6) [55]. Thanks to Eqs. (33, 34) the 
practical problem of regularizing a{i) is shifted to the problem of regularizing (V'V')vac in Eq. (32) or The 
latter is given in the effective theory (6) by the logarithmically UV-divergent expression [55] 



J TT 



(27r)4 (p| + M2)2 



(35) 

reg 



Since the precise value of the cutoff is not known but only its order of magnitude, see Eq. (7), it is customary 
to adjust the cutoff(s) in the corresponding regularization scheme such that the experimental values of 
(V'V')vac and in (32,35) arc reproduced [44]. In this way free parameters (cutoff, current quark mass, etc.) 
are fixed in the vacuum- and meson-sector of the effective theory (6) [44]. In this sense the xQSM - i.e. the 
baryon sector of the effective theory (6) - yields parameter-free results. (In some calculations in the xQSM 
the mass M was understood as a "free parameter" and allowed to vary in the range 350 MeV < M < 450 MeV 
[28, 29, 44]. The sensitivity of the model results to these variations was typically within (10 — 30)%. Here 
we rely on notions from the instanton vacuum model and consequently take the value M = 350 MeV which 
follows from the instanton phenomenology [54, 55].) 

Thus, our final (regularized) result for the form- factor reads (with B{t) defined in Eq. (29)) 

a{t) = mlfl B{t) ■ |l + O(MVL)}, (36) 

where m.,^ and /,r denote the physical pion mass and decay constant, respectively. It should be stressed that 
the result (36) does not follow from the interpolation formula of Ref. [46] (which would require to add the 
UV-finite discrete level contribution, cf. Ref. [63]). Instead (36) has to be considered as an approximate result 
for a{t) in the xQSM, which is justified by the parametrical smallness of the instanton packing fraction. 
In the following the parametrically small 0(M^p^^)-corrections often will not be indicated. 

V. DISCUSSION AND INTERPRETATION OF THE RESULTS 

Evaluating the final expression (36) for the form factor a{t) with the soliton profile (13) for M = 350 MeV, 
m^r = 140 MeV and = 93 MeV yields the result shown in Figs, la and lb (see Appendix B for detailed 
expressions). It should be noted that the error due to using the profile (13), instead of the self-consistent 
profile which truly minimizes the soliton energy (11), is far smaller than the theoretical accuracy in Eq. (36). 

Apart from the values of a{t) at the Cheng-Dashcn point t = 2mJ and at t = 0, and their difference 
A<j = cr(2m^) — (7(0), there is another phenomenologically interesting quantity - namely the scalar mean 
square radius related to the slope of the form factor a{t) at t = as 

a{t)=a„^(l + Url)t + 0{t^)) . (37) 



We obtain the results 
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£7^N = 67.9 MeV, (38) 

a{2ml) = 82.6 MeV, (39) 

= 14.7 MeV, (40) 

(r|) = 1.00 ftn^. (41) 
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0(t) [MeV] O(t)-O(0) [MeV] 




FIG. 1: a. The form factor cr{t) evaluated in the xQSM by means of Eq. (36) (solid line) and its dipole fit (42) as functions 
of t. b. The analytical continuation of the form factor a{t) to unphysical t > 0. The dot marks the value of <T(t) at the 
Cheng-Dashen point t = 2m^. At and above the threshold t > 4m^ the form factor is not defined, see text. 



A useful parameterization of the form factor for negative t is given by (with ctttn from Eq. (38)) 

^^^)" (1-T/M|)^ ' ^l-0-55GeV^ (42) 

The dipole fit (42) approximates (7{t) to within 2% for \t\ < 0.8 GeV^, cf. Fig. la. Thus, (7{t) decreases with 
increasing \t\ more quickly than the electromagnetic form factors where the corresponding dipole mass is 
about Mgjjj ~ 0.7GeV^ in a comparable t-region. 

The form-factor (j{t), Eq. (36), is not defined at and above the threshold t > Am\. In the vicinity of the 
threshold the form factor behaves as (see Appendix B) 

a{t) = ai In I j + 02 as t ^ 4m^ (where t < Arn^), (43) 

where oi, 02 are positive c;onstants. Interestingly, a similar divergent behaviour of o"(f) for t 4m^ is also 
observed in heavy baryon chiral perturbation theory [18-20]. There this feature arises as a peculiarity of the 
non-relativistic expansion and can be avoided by considering baryon chiral perturbation theory in manifestly 
Lorentz invariant form [21]. It is not clear whether in the xQSM this unphysical feature could also be cured 
- possibly by a more careful analytical continuation of a{t) to t > 0, e.g., by making use of (subtracted) 
dispersion relations. 

Accuracy of the approximation. Before discussing the results let us estimate the size of some of the 
contributions neglected in Eq. (36). The contribution of the discrete level in (22) is (a^N)iev = 12 MeV [63] 

(cf . next paragraph) . The contribution to the continuum part of cr,rN from the second order of the gradient 

(2) 

expansion is (cr7rN)co„t — —6 MeV (cf. App. A). These are corrections of 0(15%) to (38), i.e. smaller than 
the theoretical accuracy of the approximation (36) which is 0{M^pl^) = 0(30%). This is an indication (and 
no more) that the approximation works. Its theoretical justification is anyway unquestioned due to Eq. (8). 

Comparison to previous calculations in the xQSM. In [27, 28] CTttn was studied in SU(2) in leading order 
of the large- TVc expansion. In [29] the form factor was studied in the SU(3) version of the model and including 
1/A^c corrections. In [27-29] the proper time regularization was used. Our result (38) for CTttN agrees with 
the numbers quoted in [27-29] to within 30%. 

In [29] it also was observed that a{t) is not defined for t > 4:m'^. In the region < t < 2m% our result 
for the difference a{t) — (j(0) agrees with the result of Ref. [29] to to within (10 — 20)%. In the region 
t < the agreement of the rescaled form factor, i.e. a{t)/a(Q), is even more impressive (few percent). It is 
not surprizing to observe the approximation to work differently in different f-regions. The approximation 
(36) indicates the limitations of the numerical (finite box) method used in [29]. From Eqs. (B.l, B.2) in 
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Appendix B it is clear that integrals (in coordinate space) converge more and more slowly as t approaches 

t = 4m^ (from below). In the finite box method, however, it is necessary that integrals converge quickly 
in order to be in the continuum (box size — > oo) limit. Therefore, in [29] the region t > 2mJ could not be 
explored and quantitative observations of the kind (43) were not possible. 

In Ref. [63] the twist-3 distribution function e°'{x) and - by exploring the sum rule (19) - also cr,rN were 
computed in the xQSM by means of the interpolation formula of Ref. [46]. To remind, the interpolation 
formula consists in estimating the continuum contribution essentially in the same way we did here, but to 
add also the exactly evaluated discrete level contribution, which is (f7,rN)iev = 12MeV. The total result 
Cttn = 80 MeV of Ref. [63] agrees with the result (38) obtained here to within 15 %. 

To summarize, the instanton- motivated approximation for cr{t), Eq. (36), yields results in agreement with 
earlier model calculations to within the expected accuracy. 

Comparison to experimental information. The result a{2m'^) = 83 MeV in Eq. (39) agrees well with 
the more recently extracted values of a{2rn^) = (80 — 90) MeV [6-9]. The result for cr(t) — (TttN in the 
region (— 2m^) < t < 4m^, cf. Fig. lb, agrees with the result obtained from the dispersion relation analysis 
with chiral constraints of Ref. [5] to within (2 — 5)%. This can be seen by comparing the difference A^- in 
Eqs. (4, 40) and the slope of a{t) at t = 0, namely a'{0) oc ct^n (r|) = 72MeVfm^ of [5] vs. 68MeVfm^ 
obtained here. The agreement with phenomenological information can be considered as satisfactory. 

Comparison to other a,pproaches. Our result is consistent with the values for (j^rN obtained from lattice 
QCD results [23-25] (see also Section VI). In Table 1 we compare our result also to calculations performed in 
chiral perturbation theory [18, 20, 21], the linear sigma model [30], Skyrme model [31, 32], colour-dielectric 
model [33], cloudy bag model [34, 35], Nambu Jona-Lasinio model [36, 37], pcrturbativc chiral quark model 
[38], a relativistic dynamical model based on effective hadronic (A'^, A, tt, p and a) degrees of freedom [40]. 
Also the results of the data analyses and from the xQSM-calculations [27, 29] are included. 

In several models it was observed that a sizable - if not dominant - contribution to ctttn is due to the 
pion cloud ("quark sea") as compared to the constituent quark core ("valence quarks") [33-35, 38]. In the 
xQSM the discrete level contribution corresponds to the quark core and the continuum contribution to the 
pion cloud. Here we obtain an extreme picture, where cr,rN is purely due to the pion cloud. Corrections to 
this picture are suppressed by the instanton packing fraction (and are practically of order 30%). 

Skyrmion- and non- relativistic limit. It is possible to recover from expressions of the xQSM the results 
of the non-relativistic quark model and the Skyrme model by taking appropriate (non-physical) limits [68]. 

The "Skyrmion-limit" consists in taking i?soi — * oo. Since for the soliton solution Rso\ ~ [42], the 

limit is to be understood as evaluating model quantities with, e.g., the profile in Eq. (13) which allows to 
vary [68]. In this limit the energy of the discrete level i?iev {~^) [42], such that this contribution 
is enclosed in the contour of the w-integral in (25). With increasing R^oi the contributions (o',rN)corit 

series (26) behave as (,(y-Kn)^cont ^ -^soi'^- (''■7rN)cont (o-7rN)cont this can be seen directly from the 
expressions in the Appendix. For arbitrary k one arrives at this conclusion using general scaling arguments. 

Thus (o'7rN)coLt dominates again - this time, however, justified by the unphysical i?soi oo limit. The 
expression for ct^n obtained here formaJij coincides with the expressions in Refs. [31, 32]. It should be noted 
that what is an exact result in the Skyrme model is here merely an approximation - though a well justified 
one thanks to arguments from the instanton vacuum model. The coincidence of the expressions is purely 
formal since the Skyrmion is a topological soliton. E.g., in [32] a vector-meson model was used to determine 
the i7-field. 

In the opposite limit, i^soi — > 0, one recovers results from the non-relativistic constituent quark model 
(formulated for arbitrary Nc [69]). Taking i?soi ^ in the expression (16) one obtains 

a{t) = M^GeH) , (44) 

where Gsit) (normalized to Ge{0) = 1) denotes the isoscalar electric form factor - taken also in the 

non-relativistic limit. This result follows from generalizing the discussion in Ref. [63] to the case of the 
form factor a{t). That Eq. (44) is the correct non-relativistic relation between a{t) oc {N'\'tptjj\N) and 



This is to be understood in a loose sense in the "quark model language" . Strictly speaking there is no one-to-one correspon- 
dence in the xQSM between the discrete level and continuum contributions on the one hand and the notions of "valence" 
and "sea quarks" on the other hand. The latter are well defined, e.g., in the context of parton distribution functions. The 
point is that both discrete level and continuum contribute ea<;h to valence and sea quark distributions [46, 47]. 
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fT,rN/MeV a{2ml)/MeY A<,/MeV {r%Y^'^ /bo. 



data analyses 



Koch (1982) [4] 
Gasser et al. (1991) [5] 
Kaufmann and Hite (1999) [6] 
Olsson (2000) [7] 
Pavan et al. (2002) [8] 



45 ±8 



64 ±8 
60 ±8 
88± 15 
71 ±9 
79 ±7 



15.2 ±0.4 



1.3 



theoretical approaches 



conventional heavy baryon xPT[20] 
manifestly Lorentz invariant xPT [21] 



4± 1 
14 



lattice QCD, Dong et al. (1995) [23] 49.7 ± 2.6 

lattice QCD, Lcinwcbcr et al. (2000) [24] 45-55 
lattice QCD, Leinweber et al. (2003) [25] 37tll " 73t 



6.6 ± 0.6 



0.72 ± 0.09 



SU(2) xQSM to LO Nc [27] 54 

SU(3) xQSM to NLO [29] 41 59 

linear sigma model [30] 86-89 - 

SU(2) Skyrmion [31] 49 

SU(3) Skyrmion [32] 60 

chiral colour dielectric soliton [33] 37.8 - 

cloudy bag model [34] (37-47) ±9 

Nambu-Jonar-Lasinio [36] 50 - 

confined Nambu-,Jona-Lasinio [37] 60 - 

perturbative chiral quark model [38] 45 ± 5 - 

effective model of hadrons [40] - 74 



18 



1.2 



1.1 



This work (accuracy ~ 30%) 



67.9 



82.6 



14.7 



1.00 



TABLE I: Comparison of the results for iT,rN, 0"(2m^), Ao- = a(2rn^) — CTttN and (r^) from data analyses and theoretical 
approaches. The Ust is far from complete, only some of the more recent results are shown. For reviews on early approaches see 
[!]• 



GE{t) oc {N'\il;^tp\N) follows from considering that in this limit the nucleon wave-function has no lower 
(antiquark-) component such that tp^tjj = -tjjijj practically holds, and that the current quark mass m has to 
be understood as a constituent quark mass equal to Mn/3. The non-relativistic relation (T^n = strongly 
overestimates the phenomenological value of Ut^n in Eq. (5) . Still, this result is theoretically consistent and, 
e.g., correctly implies a vanishing strangeness content y in the nucleon [62]. 

Interpretation of CTttn- In the effective theory (6) the pion- nucleon sigma-term and the quark vacuum 
condensate are proportional to each other 

cr^N fX -?7l(V'?A)vac , (45) 

up to parametrically small C'(M^p^^)-corrections. From a physical point of view this observation based on 
Eq. (33) is not surprising since (V'V')vac is the sigma-term of the vacuum (per unit volume, up to the explicit 
factor of m). Thus in the xQSM the picture emerges that ct^n (measure of chiral symmetry breaking in 
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the nucleon) is directly proportional to (V'V')vac (measure of spontaneous chiral symmetry breaking) and m 

(measure of explicit chiral symmetry breaking). The relation (45) is also known from the Skyrme model 
[31, 70] which is not surprizing, see the remarks in the previous paragraph. 

The proportionality factor in Eq. (45), the function B(t) at t = as defined in (29), has the dimension of 
volume. Following the temptation to define "an effective volume of the nucleon" as VeS = B{0), the relation 
between CTttN and {iptp)va,c can be written as 

cr^N = -m ('0'0)vac V^s , Kff = B{0). (46) 

Numerically we find an efi^ective volume which - taking the nucleon to be a rigid sphere - would correspond 
to an eff'ective nucleon radius of 0.9 fm (for m„ = 140 MeV). This value should not be taken too seriously. 
But it yields the correct order of magnitude for the phenomenological size of the nucleon. Vgfi has a well 
defined chiral limit, see below. The concept of such an effective volume is useful, e.g., in the context of the 
"partial restoration of chiral symmetry in nuclear matter". The Ves = a^^N /{'m'^f^) in Eq. (46) is just the 
inverse of the chiral nucleon density introduced in Ref. [10]. 

Thus, the value of a^^ is large because (— m(V'V')vac) is sizable and because the nucleon is a large extended 
object. Corrections to this picture are suppressed by the instanton packing fraction. 

The strangeness content of the nucleon. Our result ctttN = 68 MeV implies a strangeness content in the 
nucleon of y ^ 0.4, cf. Section II, as it is also inferred from the more recent data analyses [6-9]. In view 
of such a large value for the strangeness content some authors even suggest to reconsider the "standard 
interpretation" relating f7,rN and y [8] . (However, as argued in [13] - despite the large strangeness content - 
the total contribution of the strange degree of freedom to the nucleon mass is small.) 

The calculation presented in this work definitely confirms - within its accuracy - that fT,rN is large, and 
suggests an "explanation" why ("because the nucleon is a large extended object", see previous paragraph). 
Being a calculation in a SU(2) model, however, it unfortunately cannot provide any insights into the puzzle 
why the strangeness content y appears to be so large. In this context it is worthwhile mentioning that in 
Ref. [29], where ct^n was studied in the SU(3) version of the xQSM and the matrix element {N\'tps4's\N) 
directly evaluated, the strangeness content was determined to be ~ 0.3. 

The numerical result of Ref. [29] does not provide any physical intuition why this value is so large. It 
would be interesting to see whether the present approach could be extended to the SU(3) sector and provide 
any insigths in this respect. Such a study, however, goes beyond the scope of this work. 



VI. PION MASS DEPENDENCE OF £t^n AND Mn 



The result in Eq. (36) allows us to study the form factor a{t) as a function of mir, i.e. (T(i,TO^). Of 
particular interest is hereby the behaviour in the chiral limit which will be investigated first. Maybe more 
interesting from a phenomenological point of view is the mTr-dependence of the nucleon mass Mn which can be 
deduced from 0-^^(771^) by means of the Feynman-Hellmann theorem. The understanding of the correlation 
between Mn and is presently of importance for extrapolations of lattice data - which correspond to pion 
masses of typically m,r > 500 MeV - to the physical value of the pion mass. 

Leading non-analytic contributions. Expanding the result in Eq. (36) around the point = we obtain 
for c7^N(?Ti7r), ^^('TiTr) = o'(2m^,m,r) — 0'7rN('Ti7r) and the derivative a'{0,m^) the following results 

27 

CT^N(m^) = aoml- — ^ ml + ... , (47) 



647r /2 - 
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A.(m,) = —2^mi + ... , (48) 
cT'{0,m^) = -^^m^ + ... , (49) 

where the constant ao = (3\/67r5'y^)/(32/^) and the dots denote higher order terms in the limit m,r — > 

(see Appendix B). Since tt?,^ oc to, odd powers of to^t in Eqs. (47, 48, 49) arc the respectively leading non- 
analytic contributions in the current quark mass to. Leading non-analytic contributions are of particular 
interest, because they arc model-independent. The leading non-analytic contributions in Eqs. (47, 48, 48) are 
exactly three times larger than those obtained in chiral perturbation theory. The discrepancy is explained by 
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FIG. 2: a. The ratio <T,rN(n^7r)/)n^ as function of (with mn in units of the physical pion mass and a-^tiin^Tr) in in units 
of its physical (model) value cr,rN = 68MeV, i.e. the marked point corresponds to the physical situation). The solid line is the 
full model result, the dashed line is the chiral expansion in the model up to the leading non-analytic contribution, cf. Eq. (47). 
b. MN(mx) as function of (in units of the physical pion mass). Solid line is the full result, Eq. (52). Dashed line is the 
(large- A^c) chiral expansion of Mjq{mn) up to the leading non-analytic contribution, i.e. Eq. (47) integrated by means of (51). 

recalling that here we work in the large- A^^ limit and that the limits ^ oo and m„ ^ do not commute 
[17, 71]. In the large-A^,, counting Ma — = ^i^c^) while = O(N^) such that it is appropriate to 
consider first the large- A'',, limit, which is done here (though Ma — Mjv > rn-^ in nature would suggest the 
opposite). M^{m^) has a branching point at {mq oc) = (Ma — Mn)^ which in the strict large-iVg limit 
contributes to the non-analytic behaviour of Mj^{mTr) at = 0. As a consequence the A contributes as 
intermediate state in chiral loops to leading non-analytic contributions, and its contribution is twice (for 
isoscalar quantities) that of intermediate nucleon states in the strict large Nc limit [72]. 

The result in Eq. (49) means that the scalar isoscalar mean square radius diverges in the chiral limit 

1 /2 

<^^^=(w) + ^(regular terms)}. (50) 

Another example of a square radius which diverges (however, as Inm^) in the chiral limit is the electric 
isovector charge mean square radius. (Also this feature is observed in the xQSM [44]). 

Interestingly, the correct leading non-analytic contributions follow here from the structure of the soliton, 
and not from chiral loops as in chiral perturbation theory. In Ref. [72] it was shown that the leading non- 
analytic contribution to ctttN in Eq. (47) is a general result in a large class of chiral soliton models of the 
nucleon. This result is reproduced here in the xQSM because the analytic profile (13) correctly describes 
the long distance behaviour of the chiral pion field [72] . 

Full dependence on m^r. Since (JTr'Ni^Tr) vanishes in the chiral limit as we consider the ratio 
(T^N(TO7r)/TOj- Figure 2a shows aT^j^ (niTr) / as function of mj. It is preferable to plot on the x- 
axis since the chiral limit corresponds to current quark mass m oc — » 0. The solid line shows the 
full result and the dashed line shows the expansion of (TjrN(^"7r)/"i^ i^- the model up to the leading non- 
analytic contribution, Eq. (47). Clearly, only for rather small pion masses ttItt <C 140 MeV the full result for 
o'7rN(™7r)/TO^ is reasonably approximated by its chiral expansion up to the leading non-analytic contribution 
in Eq. (47). 

Figure 2 simultaneously shows the "effective volume" of the nucleon V^fj (m^) as defined in (46) in units 
of its "physical value" (corresponding to about Sfm""*). In the chiral limit Vcs{m,-^) grows by more than 
50% compared to its physical value. But it remains finite which means that this quantity is indeed a useful 
measure of the "nucleon size" in the chiral limit - in contrast to, e.g., the scalar mean square radius (r"|), 
cf. Eq. (50). 



13 



MN(m^)/GeV M^KyGeV 




2 mJ/GeV^ 



2 mJ/GeV^ 



FIG. 3: The nucleon mass Mn as a function of according to Eq. (52) with Mn(0) = 860 MeV in comparison to (a) lattice 
QCD results from [74], and (b) from [75]. (The thick dashed line is actually the fit of Ref. [24] to the lattice results [75].) The 
meaning of the thin dotted line is explained in the text. In both figures the arrow denotes the physical point. 



Nucleon mass as function of uiTr- The Feynman-Hellmann theorem (18), reformulated by means of the 
Gell-lVIami-Oakes-Renner relation in Eq. (34) as 



cr„N(m^) = m: 



2 5Mn(to^) 



dml 



provides a mean to study the nucleon mass as function of m. 

MN(m„) = Mn(0) + / dfi 
Jo 



2 (^nNil^) 

7^ 



(51) 



(52) 



where Mn(0) is an integration constant to be identified with the value of the nucleon mass in the chiral 
limit. 

For the amount the nucleon mass is shifted in the chiral limit with respect to its physical value we obtain 
from (52) the result 



MN(140MeV) - Mi 



m = / 

Jo 



d/x^ '^'rN(M) 



= 80 MeV 



(53) 



m„=140MeV 



modulo corrections which are parametrically small in the instanton packing fraction. Within this accuracy 
the mass of the nucleon in the chiral limit is Mn(0) = 860 McV (taking the physical mass as 940 McV), which 
is in the range of the values considered in chiral perturbation theory. Since we work here in the large- A^c 
limit one obtains the same result for the mass shift of the A in the (large- A^c) chiral limit. 

Fixing the integration constant 7\ifN(0) in (52) to 860 MeV we obtain for MN(m7r) the result shown in 
Fig. 2b, where the solid line shows the full result from Eq. (52) and the dashed line shows the chiral 
expansion of MN(mTr) up to the leading non-analj^tic contribution, Eqs. (47, 52). Similarly to the case of 
a„^{m'^) the chiral expansion up to the leading non-analytic contribution reasonably approximates the full 
result for Mj^slmj^) only for below the physical pion mass. 

Several comments are in order. First, in the xQSM the nucleon mass is given by the minimum of the 
soliton energy (11) with respect to variations of the chiral field U. The accurate (and numerically involved) 
procedure to study the exact dependence of Mn on the pion mass m-rr in the xQSM would consist in 
considering the soliton energy (11) as a fimction of the current quark mass m in Eqs. (9, 10) and deducing 
the respective value of the pion mass from the Gell-Mann-Oakes-Renner relation (34), or the Yukawa- like 
decay of the self-consistent soliton profile oc exp{—m„r)/r at large r. Here, however, we can be confident to 
describe correctly the variation of Mn with m„ to within the accuracy of CTttN in Eq. (36). 
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Next, since Eq. (52) can be used to describe merely the variation of Mn with it is clear that we are 

free to choose the integration constant Mn(0) = 860 McV and that there would be no point in taking the 
precise model value for Mn(0) which, by the way, would have the drawback of introducing a regularization 
scheme dependence into the result. It should be noted that in the xQSM the baryon masses tend to be 
overestimated by about 20% due to spurious contributions of the center of mass motion of the soliton [73] . 

Finally, one has to consider that the description of M^{m'^) by means of Eq. (52) can be considered as 
well justified for rn^p^^ <C 1. This relation is analog to Eq. (8) and means that only light pions - light with 
respect to the scale - make sense as effective degrees of freedom in the low energy theory (6). 

Though, of course, it is not their actual goal, lattice QCD simulations allow to measure M-^^{m.^). The 
result for MN(m^) from Eq. (52) with M^{Q) = 860 McV is compared to results from lattice QCD simulations 
from Refs. [74, 75] in Figs. 3a and 3b. (More precisely, what is plotted in Fig. 3b as a thick dashed line is 
the parameterization of the lattice data [75] reported in Ref. [25].) 

Keeping in mind the above-mentioned reservations we observe in Fig. 3 a good agreement with the results 
for the nucleon mass from lattice QCD simulations reported in Refs. [74, 75] . A similarly good agreement is 
observed with the results reported by other lattice groups [76, 77] . 

Chiral extrapolation of lattice-QCD results. Eventually one is interested in the chiral extrapolation of 
lattice data from the currently available region 500 MeV < < several GeV to the physical value of the 
pion mass. The lattice data on M-^{m-^) arc usually fitted to ansatzc of the kind M^{mT^) = a + bra^ or 
MN(m7r) = a + hrn^ + cm^ inspired by the chiral expansion of the nucleon mass, often observing that both 
ansatze are equally acceptable [77]. However, the approximation of M-i^{mT^) by its chiral expansion up to 
the leading non-analytic term O(m^) makes sense only at small values of rriTr below the physical pion mass 

- as was emphasized in [24] (and can also be seen here in Fig. 2b). 

An interesting extrapolation method was introduced in [78] where it was suggested to regularize chiral pion 
loops (which simulate the pion cloud effect) by introducing appropriate regulators (i.e. form factors which 
simulate the extended structure of the nucleon). This approach not only incorporates the correct chiral 
behaviour, but also reproduces results of the heavy quark effective theory in the limit m cx oo. This 

gives a certain legitimation that also the intermediate rriTr-region - as explored in lattice QCD calculations - 
is reasonably described in this approach. The model dependence of this approach was studied in [24, 25] by 
using different regulators. The thick dashed curve in Fig. 3b is the fit of Ref. [25] to the lattice data [75] in 
the range 500 MeV < m^r < 1 GeV. (The thickness of the curve is comparable to the statistical error of the 
very accurate lattice data of Ref. [75].) The extrapolation from this region to the physical pion mass yields 

- within the statistical accuracy of the lattice data and depending on the regulator values for the nucleon 
mass which cover the region (782^^22 • • • 94812^6) M!eV. This demonstrates the sensitivity to details of the 
extrapolation of even very accurate lattice data. 

The approach of Refs. [24, 25, 78] basically corresponds to chiral perturbation theory in different regu- 
larization schemes. From the point of view of field theory, however, it is unsatisfactory to observe a strong 
scheme dependence [79] . Considering the complexity of the problem it is important to have further (model 
independent) constraints. Recently it was reported that chiral perturbation theory is able to describe reliably 
MN(m^) up to < 600 MeV using refined regularization techniques [79]. Presently most of the lattice 
data is beyond that limit, however, a first matching of chiral perturbation theory and lattice results seems 
to be possible [79]. 

The results reported in this work also could provide useful - though certainly not model-independent 

- insights into this issue. Indeed, by reasonably fixing the integration constant Mn(0), an agreement with 
lattice data to within the accuracy of the approach is obtained, cf. Figs. 3a and 3b. This observation suggests 
that it would be worthwhile attempting to fit lattice data with an ansatz like, e.g., 



The ansatz (54) corresponds to Eq. (52) where Mn(0) and the soliton size i?soi (cf. Eq. (13)) are allowed 
to be free parameters, Mn(O)''* and Rs^\, to be fitted to lattice data. The ansatz (54) is of more general 
character: It is not only inspired by the xQSM (and the Skyrme model, cf. previous Section V) but actually 
is based on the large Nc description of the nucleon as a chiral soliton of the pion field [43, 60]. 

The ansatz (54) means that the pion cloud effect - which is responsible for the growth of the nucleon mass 
with increasing m,r - is modelled by the structure of the soliton. This is conceptually an independent ansatz 



2 




(54) 
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to consider pion mass effects than the method of Refs. [24, 25, 78], and it has moreover the advantage of 

being regularization scheme independent. 

In order to illustrate that the ansatz (54) is reasonable note that Mn(O)''' = 780 MeV and R^^^ = 0.93i?soi 
fit the lattice data [75] within their statistical accuracy. The ansatz (54) with these parameters - shown in 
Fig. 3b as the thin dotted line in comparison to the fit obtained in [25] - yields for the physical nucleon mass 
about 850 MeV which is in the range of the values reported in Ref . [25] . A careful study whether this is the 
best fit and an estimate of its statistical and systematic errors go both beyond the scope of this work. The 
main systematic error in the ansatz (54) is due to the treatment of the A as a mass generated state to the 
nucleon in the large Nc limit and could be estimated following Ref. [72]. 

VII. SUMMARY AND CONCLUSIONS 

The sigma-term form factor a{t) of the nucleon was studied in the limit of a large number of colours in 
the framework of the chiral quark-soliton model (xQSM). The consistency of this field theoretical approach 
is illustrated by the fact that the pion-nucleon sigma-term ctttn can consistently be computed in the xQSM 
in three different ways. Apart from continuing analytically the form-factor to t = one also can make use 
of the Feynman-Hellmann theorem [56] and the sum rule for the twist-3 distribution function e"(a;) [61]. 

The model expression for a{t) was evaluated in an approximation justified by arguments from the instanton 
model of the QCD vacuum from which the xQSM was derived. The approximation is therefore theoretically 
well controlled and justified. A virtue of the approximation is that it allows to regularize i7{t), which is 
quadratically UV-divergent in the xQSM, in a regularization scheme independent way. The theoretical 
accuracy of the results is governed by the smallness of the parameter characterizing the diluteness of the 
instanton medium and is practically of 0(30%). Results from previous exact calculations in the xQSM are 
reproduced to within this accuracy. There are no adjustable parameters in the approach. 

For the form factor at the Cheng-Dashen point the value <j{2m%) = 83 MeV is found - in good agreement 
with more recent analyses of pion-nucleon scattering data [6-9]. In the region —2m^ < t < Arn^ the result 
for the form factor agrees to within few percent with the shape for a{t) obtained in [5] on the basis of a 
dispersion relation analysis. In particular cr(2m^) — cr(0) = 14.7 MeV is obtained which is close to the value 
15.2 MeV of Ref. [5]. Finally, for the pion-nucleon sigma-term the value a,rN = 68 MeV is found. 

An advantage of the present approximate result for (j{t) is its simple structure which allows, e.g., to 
study the form factor in the chiral limit and to derive the leading non-analytic (in the current quark mass) 
contributions in the model. It was shown that the expressions from the xQSM contain the correct leading 
non-analytic contributions in the large- A^c limit [72] . 

The results presented in this work fully confirm the observation made in other chiral models that the 
dominant contribution to cr,rN is due to the pion cloud. Indeed, here the extreme situation emerges that 
(TttN is solely due to the pion cloud - modulo small corrections suppressed by the instanton packing fraction. 
Moreover it is found that (TttN is proportional to the quark vacuum condensate, i.e. the order parameter 
characterizing spontanous chiral symmetry breaking, and to the current quark mass, i.e. the parameter 
characterizing explicit chiral symmetry breaking - a relation previously known only from the Skyrme model. 
The proportionality factor can be interpreted as an effective volume of the nucleon, which yields the physically 
appealing explanation that ct^n is so sizeable because the nucleon is a large extended object. 

The large value of CTttN implies according to the standard interpretation a surprizingiy sizeable strangeness 
content of the nucleon of about y = 2{N\'tl>siJsW)/{N\tl>u'fpu + i'dtpd\N) ~ 0.4. The SU(2) calculation 
presented here unfortunately cannot shed any light on this puzzle. Calculations in the SU(3) version of 
the xQSM, however, naturally accomodate a large strangeness content [29]. Worthwhile mentioning in this 
context is the conclusion [13] that despite the large value of y the total contribution of strange quarks to the 
nucleon mass is small. 

The dependence of the nucleon mass Mn on the pion mass was obtained from (T7rN(™7r) by exploring the 
Feynman-Hellmann theorem. In the chiral limit the nucleon mass was found to be reduced by 80 MeV with 
respect to its physical value. This means that the mass of the nucleon in the chiral limit is about 860 MeV 
which is in the range of the values considered in chiral perturbation theory. To within the accuracy of the 
approach the functional dependence of the nucleon mass on the pion mass agrees with lattice QCD results 
up to the largest available lattice values of m,r- This observation can be used to inspire ansatze for the chiral 
extrapolation of lattice QCD results. 

Despite the considerable progress of the rigorous approaches to QCD - lattice QCD and the effective 
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field theory method of chiral perturbation theory - the description of the pion-nucleon sigma-term still is a 

demanding issue. The direct calculation of ct^n on the lattice meets the problem that the corresponding op- 
erator is not renormalization scale invariant [23] . This complication is avoided by the method of determining 
ctttN from the pion mass dependence of the nucleon mass [74-77]. In either case, however, the conceptual 
problem appears of how to extrapolate lattice data from the presently available m^- > 500 MeV down to the 
physical value of the pion mass. A model-independent guideline can, in principle, be delivered by the chiral 
perturbation theory where, however, the physical value of ct^n cannot directly be computed since it serves 
to absorbe counter terms and is renormalized anew in each order of the chiral expansion. 

Recently a first promising matching of chiral perturbation theory to the lowest presently available pion 
masses in lattice QCD was reported [79]. However, it will take some time to demonstrate and establish the 
convergence of the chiral expansion up to values of TOtt ~ 500 MeV and higher. Until then less rigorous 
but phenomenologically and theoretically well motivated approaches may have a chance to serve as helpful 
guidelines for the chiral extrapolation of lattice data. One such approach, suggested in Ref . [78] , consists in 
exploring physically motivated regularization technics to simulate pion cloud effects. In order to investigate 
the model dependence it would be desirable to consider also other conceptually different approaches. 

The observations made in this work inspire an alternative, regularization scheme independent, extrapola- 
tion ansatz in which the pion cloud is modelled on the basis of the soliton picture of the nucleon. In view of 
the success of this picture in describing baryon properties siich ansatze appear to be promising. To provide 
a theoretically well controlled extrapolation ansatz it will, however, be necessary to introduce systematically 
1/A^c-corrections which is subject to current investigations. 

It would be interesting to study also other observables of the nucleon in this way, using arguments from the 
instanton vacuum model - which not only largely simplifies the calculations but is crucial, e.g., to allow the 
derivation of leading non-analytic contributions in the xQSM. However, it is questionable whether hereby 
one could obtain also a phenomenologically satisfactory description of observables other than quadratically 
UV-divergent, and the only observable of such kind as far as the author is aware of is the scalar form factor. 



Acknowledgments 

The author thanks S. Boffi, K. Goeke and M. V. Polyakov for fruitful discussions. This work has partly 
been performed under the contract HPRN-CT-2000-00130 of the European Commission. 



APPENDIX A: GRADIENT EXPANSION 

By using + iH) = (w - iH)/{uj^ + H^) and = -p^ + + iM-/''{S7''U^^){x) in Eq. (25) we arrive 
at the series (26) with the ^(t)^^)^^. given by 



(fe) 

cont 



m Nr 



duj 

TT- Sp 
C Am 



jo{sr=i\^\)^\iv-iH) 



W2 + p2 + M2 



X ( -iM7*=(V*=C/^=)(x) 



a;2 + p2 



M2 



(A.1) 



In order to evaluate the zeroth order contribution o^{t)'f2^^ we saturate the functional trace by the complete 
set of eigenfunctions of the free momentum operator, i.e. Sp[. . .] = / ^^^^i (p|tr . . . |p) where tr denotes the 
trace over flavour- and Dirac-indices. After taking the Dirac-trace we obtain 



<^(i)cont = -mN^AM j d^x |x|) trp 



11 / 



c 



d3p 



277 J (27r)3 + p2 + Af2 



(A.2) 



The "—1" under the flavour-trace (trp) originates from the vacuum subtraction term. The w-integration 
over the contour C in (28) can replaced by an integration over the real w-axis, and we substitute uj p% to 
use the convenient Euclidean space notation. Using tv-pU = tipU^ we obtain the result in Eq. (27). 
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Similarly we obtain for the k = 1 term in (A.l) 

which is zero due to rotational invariance. 

The k = 2 term in (A.l) yields after a somehow lengthy calculation 

^(*)cont = - / d'x io(^/=^ |x|) trp U^^) (V'=C/)(x) (V^C/t)(x) . (A.4) 

For our purposes it is sufHcient to observe that cr{t)'^^^ is finite, e.g., (cwNjcont/'^ — — (3-\/2/16)A^c-^-Rsoi- 
(The equality holds for = in the profile (13).) Higher orders k >3 are also finite, i.e. 

o-(Ocont = UV-finite , k>2. (A.5) 



APPENDIX B: <T(t) AND THE CHIRAL LIMIT 



Taking the trace over flavour- indices in the expression for B(t) in (29) and inserting the result into the 
expression (36) for a{t) we obtain for the soliton profile in (13) 



oo 

a{t) = mlfl A-K J dr jo{ry/^t) (l - cos P{r)j , 



1 — cos P{r) 



2Rj^i{l + m^r)'^ cxp(-2m„r) 



(B.l) 



The form-factor can be continued analytically to the phenomenologically interesting region f > as 

f sin(rv^) 



sinh {r^/i) 



t < 0, 
t = 0, 
t > 0. 



(B.2) 



The form factor (B.l) is undefined for t > 4m^. For positive values of t the integral over r in (B.l) converges 
only for t < 4mJ because only then at large r the decay of r^[l — cos P(r)] a exp(— 2m7rr) can compensate 
the rise of joii^V—i) oc exp(r-\/<)/r. 

Chiral limit. The pion-nucleon sigma-term crn-N(?Tin-) as function of m^- is given by 



'dx (1 -|- ax)"^ exp{—2ax) 



/dx X [1 



The zeroth order in the Taylor of the function J(a) around a = is 

oo 



(B.3) 



2^2 



(B.4) 



The linear term in the Taylor expansion of I{a) requires careful treatment. Consider a ^ and introduce an 
upper limit D ^ JR^oi in the space integral in Eq. (B.l). This simulates the "realistic situation" of computing 
(TttN in a finite volume (as it happens in many model calculations and in lattice QCD). Then 



dlja) 
da 



= -2a 



dx 2:^(1 -|- ax) exp(2aa;) 
(x4exp(2aa;) + {1 + ax)"^)"^ 



Dm„ 



= -2 



dz z^{l + z) exp(2z) 
(z4exp(2z)-Fa4(l-Fz)2)2 



(B.5) 
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where we substituted z ■ 
first, and only then m.^ 



ax in the second step and used a = m^-Rgoi- Taking the continuum limit D 
* (apperantly these limits do not commute) we arrive at 



oo 



dlja) 
da 



oo 

'■I 



2 dz {1 + z) ex.p{-2z) 



(B.6) 



Thus, we obtain 



TT 3 

I (a) = — p — 7^ a + ©(higher orders in a) 
2v 2 2 



(B.7) 



Inserting the result (B.7) into Eq. (B.3) and making use of the relation (14) in order to eliminate i?soi in 
favour of the pion decay constant /tt and nucleon axial coupling constant qa yields the result in Eq. (47). The 
results in Eqs. (48, 49) follow in an analogous way. It is not possible to proceed with the Taylor expansion 
of I{a) to still higher orders of a. This is possibly due to the fact that the next contribution in the chiral 
expansion of I{a) is oc a^lna, as it is in chiral perturbation theory. 
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